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ABSTRACT 

The stability properties of relativistic stars against gravitational collapse to black hole is a 
classical problem in general relativity. A sufficient criterion for secular instability was es- 
tablished by Friedman, Ipser and Sorkin (1988), who proved that a sequence of uniformly 
rotating barotropic stars is secularly unstable on one side of a turning point and then argued 
that a stronger result should hold: that the sequence should be stable on the opposite side, with 
the turning point marking the onset of secular instability. We show here that this expectation 
is not met. By computing in full general relativity the F-mode frequency for a large number 
of rotating stars, we show that the neutral-stability point, i.e., where the frequency becomes 
zero, differs from the turning point for rotating stars. Using numerical simulations we validate 
that the new criterion can be used to assess the dynamical stability of relativistic rotating stars. 

Key words: relativistic processes - methods: numerical - stars: neutron - stars: oscillations 
- stars: rotation - black hole physics 



■ 1 INTRODUCTION 

, The stability of a relativistic star against collapse to black hole 
• is one of the most important predictions of general relativity. 
' Whil e this problem is re asonably well understood for nonrotating 
, stars ( lMisneretalll973l) . this is not the case for rotating stars and 
' is particularly obscure when the stars are rapidly rotating. A mile- 
I stone in this landscape is the criterion for secular stability proposed 
. by Friedman, Ipser and Sorkin (1988), who proved that a sequence 
' of uniformly rotating barotropic stars is secularly unstable on one 
, side of a turning point (an extremum of mass along a sequence of 
i: constant angular momentum, or an extremum of angular momen- 

■ tum along a sequence of constant rest-mass). They then argued, 
based on an expectation that viscosity leads to uniform rotation, 
that the turning point should identify the onset of secular instabil- 
ity. While for nonrotating star the turning point coincides with the 
secular-instability point (and with the dynamical-instability point 
for a barotropic star if the perturbation satisfies the same equation 
of state of the equilibrium model), for rotating stars it is only a 
sufficient condition for a secular instability. Lacking other guides, 
the turning point is ro utinely used to find a dynami cal instability in 
simulations jBaiotti et al. 2005; Rad ice et al]|2010h . 

Our understanding of the dynamical instability of relativis- 
tic stars in uniform rotation can be improved by determining the 
neutral-stability line, that is the set of stellar models whose fre- 
quency of the fundamental mode of quasi-radial oscillation (F 
mode) is vanishingly small. While this problem is challenging from 
a perturbative point of view, especially when the rate of rotation be- 
comes high, it can be tackled through numerical calculations. We 
have therefore simulated in full general relativity 54 stellar mod- 
els and calculated accurately the corresponding F-mode frequency 



via a novel analysis of the power spectral density (PSD) of the 
central rest-mass density. This new approach has been validated 
through a comparison with all the available data, showing excellent 
agreement and, most importantly, a much smaller variance. By con- 
struction, in fact, simulations cannot evolve models at (or near) the 
neutral-stability line, but the accuracy of our _F-mode frequencies 
and their smooth dependence on the central rest-mass density and 
dimensionless rotation rate, have allowed us to produce an analytic 
fit of the data and deduce from this the neutral-stability line. We 
find in this way that it coincides with the turning point for spherical 
stars, but not for rotating stars, with the difference increasing with 
the angular momentum. Although somewhat suiprising, this dif- 
ference is not in contrast with the predictions of the turning-point 
criterion, since the latter is only a sufficient condition for secular 
instability and not a necessary condition for secular and dynamical 
instability. Hence, a stellar model which is stable according to the 
turning-point criterion, can be nevertheless dynamical unstable. 

To test the new stability line and validate whether it can be 
used to mark the threshold for dynamical stability, we have evolved 
stellar models whose properties fall in a small region near the two 
stability lines. Special attention has been paid to stellar models that 
are predicted to be stable by the turning-point criterion but unstable 
by the neutral-stability line. Because these model indeed collapse 
to black holes, we conclude that the neutral-stability line can be 
used effectively to mark the boundary to dynamical instability. 

The organization of the paper is as follows. Section|2]describes 
the numerical setup and initial data, while Sect.|3]presents our ap- 
proach to extract the eigenfrequency and offers comparisons with 
previous work. Section|4]collects our results and a comparison be- 
tween the two stability criteria, leaving the conclusions to Sect. [5] 
Unless stated differently, we use units in which c — G = Mq = 1. 
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2 NUMERICAL SETUP AND INITIAL DATA 

All of our calculations have been performed in full general 
re lativity (GR) using th e Whisky 2D code described in detail 
in lKellerman et aljj2008l) . This is a 2-dim ensional (2D) code b ased 
on the 3-dimensional (3D) Whisky code l lBaiotti et al.ll200^ . and 
exploitin g the condition of ax isymmetry through the "cartoon" 
method JAlcubierre et al. II2OOI]). In essence, the evolution of the 
spacetime is obtained using the 2D version of Ccatie, a finite- 
differencing code providing the solution o f a conformal trace - 
less formulation of the Einstein equations jPollnev et alj l2007h . 
while the equations of relativistic hydrodynamics are solved a flux- 
conse rvative formulation of the equations, as first discussed in de- 
tail in lBaiotti et al. The Whisky2D code implements a va- 
riety of approximate Riemann solvers and severa l recon struction 
methods and, as discussed in iGiacomazzo et al. I l l2009l) , the use 
of reconstruction schemes of order high enough is fundamental 
for an accurate evolution. In particular, the results presented here 
have been com puted using the piecewise-pa rabolic reconstruction 
method(PPM) jColella & Woodwiu-dl 1 19841) . the HLLE approxi- 
mate Riemann solver 1 Harten et al ]ll983h , and a 3rd-order Runge- 
Kutta method for the time evolution. 

The initial equilibrium stellar m odels are built using the rns 
code dStergioulas & Friedmanll9^5h as isentropic, uniformly rotat- 
ing relativistic perfect-fluid polytropes with equation of state 



p = hp 



e = p- 



r - 1 



(1) 



where p is the pressure, p the rest-mass density, K the polytropic 
constant, F the polytropic exponent, and e the energy density. Al- 
though all the results can be rescaled for any choice of K and F, 
we have here set K = 100 and F = 2, which yield stars with max- 
imum gravitational mass is M = 1.64 Af© for a nonrotating star 
and M = 1.88 A/q for a uniformly rotating one. The rns code 
provides an equilibrium solution in spherical polar coordinates af- 
ter specifying for each stellar model a central density pc and a equa- 
torial and polar (coordinate) radii in a ratio rp/re- Once this solu- 
tion is found, it is mapped to a Cartesian grid of Whisky2D and 
used as initial data for the subsequent evolution. Attention needs 
to be paid that the resolution in the calculation of the initial data 
matches well the one used in the evolution. We have verified that a 
resolution of (n^.ne) = (2001,2601) [(jir.ne) = (1001,1301)], 
with (rir, rig) the number of points of the radial and angular grids 
of the rns code, are needed for an accurate evolution in the high 
[low] -resolution setup of the Whisky2D code. Furthermore, be- 
cause we are not interested here in extracting gravitational-wave in- 
formation, we place the outer boundary at a few stellar radii and use 
a unifom grid with spacing Ax = Az = h ranging between h = 
0.04 Mq for the rapidly rotating mode ls and up to = 0.1 M p, for 
the slowly rotating ones. As done in iKellerman et"ai] ( |2008|) . we 
stagger the grid in the a::-direction of half a cell. A large number of 
tests have been carried out to verify that the results do not depend 
on the position of the o uter boundary , or on the value of the density 
in the atmosphere (see lBaiotti et alj ( |2005|) ). which we set to be 9 

orders of magnitude smaller than the cent ral one. 

As discussed by many authors jFont et al.l I2OO0I |2002| : 
iBaiotti et alJl2Q05h . the truncation error in the initial data is suf- 
ficient to trigger perturbations in the star, which will start to os- 
cillate in a number of eigenmodes. However, because we need to 
determine the eigenfrequency of the F mode, it is important that as 
much as possible of the initial perturbation energy goes into excit- 
ing that mode. For this reason we introduce an initial perturbation 
using the eigenfunction of the F mode for a nonrotating neutron 



star with the same central density, and which can be computed from 
linear perturbation theory. More specifically, denoting with V'tov 
any fluid quantity of the nonrotating model with the same central 
density and with (5V'tov ('") '^^e corresponding eigenfunction with r 
the radial coordinate in isotropic coordinate system, we as approx- 
imate equivalent eigenfunction for a rotating star in a coordinate 
system {r,6) as 5tp{r,9) = S^^^^{rR^^^ / R{9)), where R-j.^^^ 
is the radius of the nonrotating star and R{9) that of the rotating 
star, which will obviously depend on the angle ^. As a result, the 
power in the initial perturbation is mostly concentrated in the F 
mode, whose corresponding peak in the PSD of any hydrodynam- 
ical quantity is larger by at least a factor 10 than any other mode. 
As an additional validation of the procedure, we have computed the 
numerical eigenfunction for some selected models and verified that 
it matches very well the guessed one even in the case of rapidly 
rotating stars and long-term evolutions. 



3 METHODOLOGY AND ACCURACY 

As customary, we extract the f -mode frequency by perform- 
ing a discrete Fourier transform of the evolution of a representa- 
tive hydrodynamical quantity, such as the central rest-mass den- 
sity Pc and by inspecting the coiTesponding PSD. Defining as 
Fn the frequency of the largest peak in the numerical PSD, pre- 
vious studies determined the value of the F-mode frequency, 
F, by fittin g the PSD with a kn own analytic function [e.g., a 
Lorentzian, ( Kellerman et al.ll2008h ] or by taking the derivative of 
the PSD jZink et alj|2010l) . The frequency obtained in these ways 
depends sensitively on the fitting function used, on the shape of the 
PSD around Fn, and on the evolution time r. We here use a dif- 
ferent approach. Because Fn will tend to F as the evolution time 
T — >■ 00, we simply consider the evolution of Fn for increasingly 
large values of r. What we find in this way is that Fj^{t) is an os- 
cillating function around F, whose amplitude is however bounded 
by two envelopes which have a clear 1/r dependence. Fitting for 
these envelopes and extrapolating for r — > 00 we obtain a very 
accurate and possibly optimal value for F. As we will discuss in 
the following Section, this approach turns out to give an excellent 
measure of the F-mode eigenfrequency and we recommend it in all 
those studies aimed at determining eigenfrequencies of relativistic 
stars. 



3.1 Comparison with previous works 

The F-mode frequency of spherical stars can be computed 
to arbitrary precision within a linear perturbative approach 
(see l Yoshida & Eriguchil 1 I2OOII) and references therein). Hence, as 
a first validation of the accuracy of our procedure we have es- 
timated the F-mode frequency for 14 nonrotating models with 
Pc e [3.0 X 10-*, 3.0 X 10^=^]; in this range the F mode first 
grows, then reaches a maximum, and finally decreases to zero at 
the secular-instability point, around pc ~ 3.18 x 10^''. Defin- 



ing the relative error as a^, 



2 

PT 



{F)l]/{F)l^, where 



Fp.p and F^ are respectively the frequencies of the F mode from 
perturbation theory and from our simulation. The relative error is 
extremely small at low densities (e.g., (Jrci ^ 0.005 for pc ~ 
0.3— 1.5 X 10^'^) and it increases with the density (e.g., a^^i ^ 0.05 
for pc ~ 2.5 — 3.0 X 10^^), becoming of the order of about 10% at 
the edge of the secular instability. This is obviously due to the fact 
that as Fn ~ 0, numerical calculations become increasingly long 
and inaccurate. 
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Figure 1. Comparison of our i^-mode frequencies with those of 
previous wor ks in either pertu rb ation theory (PT), the CFC ap- 
prox imation jpi mmelmeier et al.l |2006|), th e Cowling approxima- 
tion jOaertig & Kokkotadl2008l : lzink et ainOlOh . or in full GR izink et alj 



We next compare our numerical estimates for the F-mode 
frequency with those made in several different approaches and 
approximations, using as reference a central rest-mass den- 
sity pc = 1.28 X 10""^, as this is the one most commonly 
used. We start our comparison by considering the case of 
n onrotating stars, for which results are available from works 



of lPimmelmeier et al.l (2006 ) in eith er the conformally-flat condi- 
tion (CFC), or of IZink eTaU ( |2010|) in full GR. This is shown in 
Fig. [T] which reports the F-mode frequency as a function of the 
dimensionless ratio j3 = T/\W\ between the rotational kinetic en- 
ergy T and the binding energy W. Note that the frequency is re- 
ported in two different scales, referring to simulations either in full 
GR/CFC (left scale) or in the Cowling approximation (right scale), 
which systematically yields larger frequencies. Although the CFC 
(blue crosses) for a nonrotating should give the same frequency in 
full GR (magenta filled triangles and red crosses) and in perturba- 
tion theory (black filled circle). Fig. [T| shows that this is not quite 
the case, although the differences are only of ~ 2%. Considerably 
larger are instead the differences with the frequencies in the Cowl- 
ing approximation, which are larger of a factor of ~ 3 (green stars 
and light-blued filled squares). Clearly, the difference between the 
results in full GR and the perturbative ones is much smaller and in- 
deed the one with our new results is the smallest among all the data 
available. We also note that our results also report the estimated 
error bars, which are much smaller than the size of the symbols. 

Considering next the comparison also for rotating stars, it is 
easy to see that our results in two dim ensions match well those in 
three dimensions of IZink et al. I( l2010h for the rotation rates avail- 
able and obviously have sma ller error bars. The very g ood match 
with the results in the CFC jPimmelmeier et ai]|200d) . with dif- 
ferences of a few perce nt only for all the values o f /3, confirms 
the conclusions drawn bv lPimmelmeier et al.l ( 1200^ , that the CFC 
is a very good approximation, at least for the dynamics of iso- 
lated stars. Figure [T|also shows that the comp arison with frequen- 
cies computed in the Cowling approximation jGaertig & Kokkota^ 




PcXlO' 



Figure 2. Square of the F-mode frequencies (blue filled circles) as a func- 
tion of Pc and 13. The dashed green area shows models above the mass- 
shedding limit, and the red solid line mai'ks the neutral stability (cf. Fig.[3j. 



l2008l : IZink et alj I2OI0I) is considerably worse. Besides an in- 
tri nsic diffe rence betwe en the two sets of data (the frequencies 
of iGaertig & Kokkota^ ilOO^ are agreement only within the er- 
ror bars of Zink et alT i 2010h ). the rate of change of the frequencies 
with /3 differs from the one found in full GR, being less rapid for the 
latter (this is not evident because the figure has two different verti- 
cal scales). This comparison shows the Cowling approximation to 
be inaccurate for all rotation rates. 

In summary, this comparison validates our approach, high- 
lighting its accuracy and smoothness when compared to alternative 
methods. This will be essential to find the neutral-stability line. 



4 RESULTS 

As mentioned above, the space of parameters is spanned by cen- 
tral rest-mass density and the angular momentum of the rotating 
models. To cover the largest possible region of parameters we have 
evolved 54 stellar models of relativistic stars with pc in the rang43 
[pmin,/9max] = [8 X 10~'',3.18 X 10"'^] and dimensionless ro- 
tation parameter /3 between zero and the mass-shedding limit for 
the corresponding sequence of constant central rest-mass density 
(/3 — 0.095 is the largest value considered). In this way we com- 
puted stellar models with masses in the range AI/Mq G [1.1, 1.9]. 

We show as filled blue circles in Fig. |2] all of the computed 
F-mode frequencies, where the squares of the F-mode frequen- 
cies (F)^ are reported as function of pc and /3. Shown as a solid 
magenta line is the analytic fitting of the frequency for nonrotating 
stars, while dashed blue lines show sequences of rotating stars hav- 
ing the same rest-mass density. All models simulated have nonzero 
F-mode frequencies and their number diminishes for (F)^ ~ 0. 
As mentioned above, this is because for these models the oscil- 
lation timescale tends to become extremely large (diverging for 
F = 0), thus becoming intractable in numerical simulations. In 
addition, models near the neutral point could also be artificially 
induced to col lapse simply by the accumulation of the truncation 
error (see also IShibatal ( 120031) ), thus preventing any reliable mea- 
sure. As a result, our analysis has been constrained to values of 
the frequencies F > 2.2 x 10"^ ~ 0.45 kHz. Fortunately, how- 
ever, the quality of the data and the smoothness in which they ap- 
pear in Fig. |2] allow us to compute an analytic fit of the function 



1 Note that pc = 1.0 X 10"^ ~ 0.62 x 10^^ g/cm^ and that pmax also 
marks the secular stability point for a nonrotating star. 
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{FY = {FY{Pc, P) and thus determine analytically the neutral- 
stability line where (F)^ = 0. 

It is convenient to use a fitting function (fflt)^(pc, P) that is 
linear in j3 and such that (F)|t(Pmax, 0) = by construction 



71 = 

5 



(2) 



(3) 



where a„, 6„ are constant coefficients, which a least-square fitting 
with the data reveals to be 

as = 6.978 x 10® , a4 = -7.757 x lO'' , as = 3.621 x 10* , 
a2 = -9.599 X 10 , ai = 1.172 x 10"^ , ao = 2.110 x 10"'^ , 
fes = -5.599 X 10^° , bi = 4.862 x 10® , 63 = -1.612 x lO*' , 
62 = 2.545 X 10^ , hi = -1.896 , 60 = 3.357 x lO""* . 

A confirmation of the accuracy of the ansatz (|3) comes from 
the very small variance of a comparison with perturbative results 
for nonrotating stars. Considering in fact over 90 stellar mod- 
els with pc € [1.0 X 10"^, 3.182 x 10"^], we obtain aat = 
\iF)l^ - {F)l^{p,,0)\ < 2 x 10-^ ~ 8 X 10-=* (kHz)2. Simi- 
larly, when comparing over the whole set of numerical data we find 
a variance that, as expected, is greater for large values of pc and 
/3 but that, overall, is aat < (Tmax ~ 1 x 10~^. Note that these 
errors are smaller or at most comparable with the numerical error 
bar, highlighting the quality of the fit. 

Using expression ([3j, it is straightforward to compute the 
neutral-stability line in a (pc, /3) plane as the one at which 
(^)flt (Pc, /3) = 0. Of course this line will be "thickened" by the 
uncertainty associated to the fit which, to be conservative, we con- 
sider to be (Jmax (We note that the thickness is much smaller for 
/? ~ but it may be larger at high /? as a result of the extrapola- 
tion.). While a neutral-stability line is already very informative in 
a (pc P) plane, its greatest impact can be appreciated in the more 
traditional (pc, M) diagram. This is shown in Fig. [3] where the two 
solid black lines refer to sequences of nonrotating (lower line) and 
mass-shedding models (upper line), respectively. Drawn as solid 
red is the neutral-stability line "thickened" by the error bar (Jmax 
(black dot-dashed lines). Finally, shown as a blue dashed line is 
the turning-point criterion for secular stability along a sequence of 
constant angular momentum J, i.e., [dM /dpc) j^^^^^^ = 0. 

Clearly, the new neutral-stability criterion does coincide with 
the turning-point criterion for nonrotating stars (cf. small inset), but 
it differs from it as the angular momentum is increased, moving to 
smaller central rest-mass densities. While unexpected, this differ- 
ence does not point to a conflict between the two criteria. This is 
because the turning-point criterion is only a sufficient condition for 
secular instability of rotating stars; stated differently, while a ro- 
tating stellar model which is at or to the right of the turning-point 
line is expected to be also secular unstable, the opposite is not true. 
Hence, the two criteria are compatible as long as the secular insta- 
bility line lies to the left (i.e., for smaller central rest-mass densi- 
ties) of the neutral-stability line. Determining the secular-stability 
line requires to consider a dissipative mechanism such as viscos- 
ity, which is however absent in our perfect-fluid description and 
difficult to introduce within a fully relativistic hyperbolic descrip- 
tion. However, because a dynamically unstable model should also 
be secularly unstable, we in fact expect the secular stability line to 
coincide or to be on the left of the neutral-stability line. In other 
words, along a J = const, sequence of stellar models we expect 
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Figure 3. Stability lines in a (pc, M) diagram. The two solid black lines 
mark sequences with either zero (lower line) or mass-shedding angular mo- 
mentum (upper Une), with the filled symbols marking the corresponding 
maximum masses. The solid red line is the neutral-stability line, "thick- 
ened" by the error bar (black dot-dashed lines). The blue dashed line is 
instead the tuniing-point criterion for secular stability. Marked with empty 
or filled circles are representative models with constant angular velocity 1 , 
02, 03, or constant initial central rest-mass density Rl, R2, R3. 



the following order with increasing rest-mass density: secular in- 
stability, dynamical instability, turning-point. 

To validate that the neutral-stability line should be used in 
place of the turning-point line to distinguish stellar models that are 
dynamically unstable from those that are instead stable, we have 
considered 6 representative models whose properties fall in a small 
region near the two stability lines. More specifically, we consider 
two different sequences having either constant angular velocity, 
i.e., models 01, 02, 03 in Fig.|4l or constant pc, i.e., models Rl, 
R2, R3. The predictions for these models are different according to 
which criterion is used for stability. In fact, while models 01, Rl 
are expected to be stable for both criteria and models 03, R3 are 
expected to be unstable for both criteria, models 02, R2 are pre- 
dicted to be stable on a dynamical timescale by the turning-point 
criterion but unstable by the neutral-stability criterion. 

To test these predictions we have evolved these configurations 
maintaining the same computational setup (but without an initial 
perturbation) and collected the coiTesponding evolution of the cen- 
tral rest-mass density in Fig. |3] As expected, models 01, Rl are 
found to be stable over about 7 ms as indicated by the central rest- 
mass density that remains constant (modulo the F-mode oscilla- 
tions), while models 03, R3 are found to collapse to black holes in 
less than 2 ms as indic ated by the exponen ti al increase of the res t- 
mass density (see also iBaiotti et al.l HoO^ : iRadice et alj told) ). 
Similarly, models 02, R2 are also found to collapse to black holes 
over a timescale which is only slightly larger than that of mod- 
els 03, R3. After validating that these results do not depend on 
the specific numerical setup used {e.g., placement of outer bound- 
aries, resolution or density in the atmosphere), we conclude that the 
neutral-stability line can indeed be used to mark the boundary of a 
dynamically unstable region. 
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Figure 4. Evolution of pc for models with constant angular velocity (upper 
panel) or constant initial central rest-mass density (lower panel). An expo- 
nential growth signals the collapse to black hole (c/ Fig.|3)- 



cial attention to those stellar models that are predicted to be sta- 
ble on a dynamical timescale by the turning-point criterion but un- 
stable by the neutral-stability line. Numerical evidence that these 
model do collapse to black holes allows us to conclude that the 
neutral-stability line can be used effectively to mark the boundary 
to dynamical instability. Besides improving our understanding of 
the stability of relativistic stars, these results show that producing 
black holes via the gravitational collapse of a neutron star is sim- 
pler than expected. Furthermore, they can serve as a guide when 
determining the neutral-stability line via perturbative techniques or 
when extending it to differentially rotating stars. 
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5 CONCLUSIONS 

The stability of rotating relativistic stars against gravitational col- 
lapse to black hole is an old problem in general relativity, impacting 
all those astrophysical problems where a neutron star may be pro- 
duced and induced to collapse as a result of mass accretion. Despite 
the importance of this problem, no analytic a criterion is known 
for a dynamical stability of rotating stars. Important progress was 
made about 20 years ago, when a cr iterion for secular stability was 
proposed by [Friedman et al.l ( Il988h . who suggested that a turning 
point along a sequence of stellar models with constant angular mo- 
mentum can be associated with the onset of a secular instability. 
Although this criterion is only a sufficient condition for the devel- 
opment of a secular instability, it has been systematically used to 
limit the region of dynamical instability in sim ulations of relativis- 
tic stars teaiotti et al.l2005l : lRadice et al.ll201(]h . 

To improve our understanding of the dynamical instability of 
relativistic stars in uniform rotation, we have computed the neutral- 
stability point for a large class of stellar models, i.e., the set of stel- 
lar models whose F-mode frequency is vanishingly small (in a non- 
rotating star this point marks the dynamical stability limit). More 
specifically, we have evolved in full general relativity 54 stellar 
models and calculated the corresponding F-mode frequency via a 
novel analysis of the PSD of the central rest-mass density. Although 
our simulations cannot evolve models near the neutral-stability line, 
the high accuracy of our estimates for the eigenfrequencies (which 
have been validated through a comparison with all the available 
data) and their regular dependence on the central rest-mass density 
and dimensionless rotation rate, have allowed us to produce an an- 
alytic fit of the data and deduce from this the neutr al-stability line. 
The la tter coincides with the turning-point line of [Friedman et al.l 
( Il988l) for nonrotating stars, but differs from it as the angular mo- 
mentum is increased, being located at smaller central rest-mass 
densities as the angular momentum is increased. This difference 
does not contradict turning-point criterion since the latter is only a 
sufficient condition for secular instability. 

To test this result we have evolved stellar models whose prop- 
erties fall in a small region near the two stability lines, paying spe- 
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